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ABSTRACT 

Ground-based 2.8-4.1 /im slit spectra of the nuclei of seven ultraluminous 
infrared galaxies (ULIRGs) that are classified optically as LINERs and have cool 
far-infrared colors are presented. All the nuclei show 3.3 polycyclic aro- 
matic hydrocarbon (PAH) emission, with equivalent widths that are systemati- 
cally lower than those in starburst galaxies. Strong 3.1 /im H 2 ice absorption, 
with optical depth greater than 0.6, is also detected in five nuclei, and 3.4 /im 
carbonaceous dust absorption is detected clearly in one of the five nuclei. It is 
quantitatively demonstrated that the large optical depths of the H 2 ice absorp- 
tion in the five sources, and the 3.4 /im absorption in one source, are incompatible 
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with a geometry in which the energy sources are spatially mixed with dust and 
molecular g clS, clS IS expected for a typical starburst, but instead require that a 
large amount of nuclear dust (including ice-covered grains) and molecular gas be 
distributed in a screen in front of the 3-4 /im continuum-emitting sources. This 
geometrical requirement can naturally be met if the energy sources are more cen- 
trally concentrated than the nuclear dust and molecular gas. The low equivalent 
widths of the PAH emission compared to starbursts and the central concentra- 
tion of the nuclear energy sources in these five ULIRGs are best explained by the 
presence of energetically important active galactic nuclei deeply buried in dust 
and molecular gas. 

Subject headings: galaxies: active — galaxies: nuclei - - infrared: galaxies - 
galaxies: individual (UGC 5101, IRAS 00188-0856, IRAS 03250+1606, IRAS 
16487+5447, IRAS 21329-2346, IRAS 23234+0946, and IRAS 23327+2913) 

1. Introduction 

Galaxies that radiate most of their extremely large, quasar-like luminosities 10 12 L Q ) 
as infrared dust emission - the ultraluminous infrared galaxies (ULIRGs; Sanders & Mirabel 
1996) - dominate the bright end of the galaxy luminosity function in the nearby universe 
(Soifer et al. 1987). They have been used, extensively, to derive information on the dust- 
obscured star-formation rate, dust content, and metallicity in the early universe (Barger, 
Cowie, & Richards 2000; Blain et al. 1999). Understanding the nature of nearby ULIRGs, 
and more particularly, determining whether they are powered by starbursts or active galactic 
nuclei (AGNs), or both, is of great importance in modern extragalactic astronomy. 

Spectroscopic observation of the thermal infrared wavelength range (3-20 /im) is cur- 
rently one of the most powerful methods of determining the energy sources of ULIRGs. At 
this wavelength range, dust extinction is lower than at shorter wavelengths (<2 /im), so 
it becomes possible to detect and measure the emission from more highly obscured energy 
sources, with smaller uncertainty in dust extinction corrections. In addition, spectral features 
in this waveband can be used to distinguish between starburst and AGN activity. Polycyclic 
aromatic hydrocarbon (PAH) emission is detected in starbursts but not in AGNs, making it 
a good indicator of starburst activity (Moorwood 1986; Roche et al. 1991), while the presence 
of AGNs hidden behind dust can be recognized through absorption features. PAH emission 
and absorption features, detected in ground-based slit spectra (e.g., Dudley 1999; Imanishi 
& Dudley 2000; Soifer et al. 2002) and in ISO spectra taken with large apertures (e.g., Tran 
et al. 2001), have been utilized to investigate the energy sources of ULIRGs. 
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In ULIRGs, there are two potentially energetically significant components. One is a 
weakly obscured (A v ^5 20 mag) 3 starburst in the extended (kpc scale) host galaxy and 
also at the nucleus. The other is a nuclear, compact (less than a few hundred pc in size), 
highly obscured energy source, powered by an AGN and/or a compact starburst. Recent 
ground-based high spatial resolution imaging observations (Surace & Sanders f999; Surace, 
Sanders, & Evans 2000) have shown that the weakly obscured (Ay < 20 mag) starbursts are 
energetically insignificant, and that the nuclear, compact, highly obscured energy sources 
dominate in ULIRGs (Sakamoto et al. 1999; Soifer et al. 2000). However, the nature of these 
nuclear sources is not easy to determine, because dust extinction can be significant even in 
the thermal infrared wavelength range. The observed PAH-to-infrared luminosity ratios in 
ULIRGs have been found to be smaller than in known starburst galaxies by a large factor 
(Fischer 2000). Assuming that the luminosities of the PAH emission are a good indicator of 
the absolute magnitudes of starbursts, the smaller ratios in ULIRGs could be due either to an 
increase in the dust extinction of the PAH emission, or a dominant contribution from AGNs 
to the infrared luminosities of ULIRGs. Determining which explanation is correct requires 
reliable estimates of the dust extinction toward the starbursts in ULIRGs. However, these 
estimates become complicated if the assumption of a "foreground screen" dust geometry is 
not valid for the starbursts (McLeod et al. 1993). 

To understand the nature of the compact, highly obscured energy sources in ULIRGs, 
the spectral shapes of the nuclear compact emission can provide more insight than just the 
absolute luminosities of PAH emission, because the equivalent widths of PAH emission are, 
by definition, robust to the amount of dust extinction. If the nuclear compact emission 
of ULIRGs displayed significantly lower equivalent widths of PAH emission than typically 
seen in starburst galaxies, and, in addition, strong absorption features were detected, it 
would suggest that energetically important dust-obscured AGNs were present. Since the 
attenuation of the extended, weakly obscured emission is much smaller than that of the 
emission from the compact, highly obscured energy sources, the observed 3-20 /xm spectra 
taken with large apertures can be strongly contaminated by the extended emission, even 
though the total infrared (8-1000 fxm) luminosities of ULIRGs are dominated by the compact 
energy sources. In this case, the signatures of absorption features toward the compact energy 
sources can be diluted. Hence spectroscopy with a narrow (less than a few arcsec) slit - 
comparable to the spatial extent of the compact component - is better suited to probe the 
nature of the compact energy sources in ULIRGs. 



3 Throughout this paper, "weakly obscured" is used to mean that the flux attenuation at 3-4 /xm is 
insignificant, which is roughly the criterion that A3_4 Mm ^ 1 mag, or Ay i5 20 mag if the Galactic extinction 
curve of A 3 -4 lim /Ay ~ 0.06 (Rieke & Lebofsky 1985; Lutz et al. 1996) is adopted. 
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Ground-based slit spectroscopy in the L-band (2.8-4.1 /im) has several important ad- 
vantages for investigating the nature of dust-enshrouded AGNs: 

1. The dust extinction curve is fairly grey (wavelength-independent) at 3-10 /iin (Lutz et 
al. 1996). 

2. The dust around an AGN has a strong temperature gradient, in that the inner (outer) 
dust has a higher (lower) temperature. The temperature of the innermost dust is de- 
termined by the dust sublimation temperature (T ~ 1000 K). In the case of blackbody 
radiation (A pca k [/im] x T [K] ~ 3000), the 3-4 /im continuum emission is dominated 
by dust with T = 800-1000 K located in the innermost regions, while the ~10 /im 
continuum emission is dominated by dust with T ~ 300 K located further out. Thus, 
the dust extinction toward the ~10 /im continuum emission regions, as estimated using 
e.g., the optical depth of the 9.7 /im silicate dust absorption feature, is much smaller 
than the dust extinction toward the AGN; the dereddened AGN luminosity will be 
underestimated if the extinction correction is made based on ~10 /im data alone. By 
contrast, the dust extinction toward the innermost, 3-4 /im continuum emission re- 
gions, as estimated using 3-4 /im data, is a better indicator of the dust extinction 
toward the AGN itself. In the case that an AGN is the only energy source, the 3-4 /im 
continuum emission predominantly originates in the innermost, 800-1000 K dust. The 
dust extinction toward the 3-4 fim continuum emission can therefore be reliably fit- 
ted with a foreground screen dust model. In a real galaxy, weakly obscured starburst 
emission may contribute to the observed nuclear 3-4 fim flux. However, as long as 
the AGN is sufficiently luminous that AGN-powered continuum emission contributes 
significantly to the observed 3-4 /im flux, evidence for an AGN can be found, pro- 
vided that the starburst contribution is subtracted properly, even if the absorption 
optical depth toward the AGN is larger than unity at 3-4 /tm (A 3 _ 4/jm > 1 mag or 
Ay > 20 mag; Rieke & Lebofsky 1985). We have previously demonstrated this in two 
ULIRGs: IRAS 08572+3915 (Ay > 130 mag; Imanishi & Dudley 2000), and UGC 
5101 (Ay > 115 mag; Imanishi et al. 2001). 

3. At 2.8-4.1 /tm, in addition to the 3.3 /im PAH emission, 3.4 fim carbonaceous dust 
absorption and 3.1 fim H2O ice absorption are present. As we will show later in this 
paper, all these features can be used in a complementary manner to constrain the 
properties of the compact energy sources. They are simultaneously observable with 
recently available spectrographs that cover the whole L-band atmospheric window (2.8- 
4.1 /tm), if the target sources are at z < 0.15. In this redshift range, a statistically 
significant number of ULIRGs have been found with IRAS. 
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4. In the near future, slit spectra of ULIRGs at 5.3-40 fim will be obtainable using the 
SIRTF IRS. By combining 3-4 fim and 5.3-40 /im slit spectra, the presence of a strong 
temperature gradient in the dust (a signature of an AGN) can be investigated through 
the comparison of the optical depths of absorption features at different wavelengths at 
3-20 fim (Dudley k Wynn- Williams 1997; Imanishi 2000). 

In the study of ULIRGs, it is particularly important to detect and quantify elusive 
AGNs that are deeply buried in dust along all sightlines (hereafter buried AGNs; Imanishi 
et al. 2001). ULIRGs that are classified optically as LINERs and have cool far-infrared 
colors (Sanders et al. 1988) have generally been taken to be starburst-powered (Downes 
& Solomon 1998; Taniguchi et al. 1999; Lutz et al. 1999). However, Imanishi et al. (2001) 
argued that LINER ULIRGs may contain buried AGNs (see also Antonucci 2001). In this 
paper 2.8-4.1 /im ground-based slit spectra of LINER ULIRGs with cool far-infrared colors 
are reported, to search for observational evidence for energetically important, buried AGNs. 
Throughout this paper, H = 75 km s _1 Mpc -1 , Qm = 0.3, and = 0.7 are adopted. 

2. Targets, Observations and Data Analysis 

The seven LINER ULIRGs in Table 1 were observed. Table 1 summarizes their infrared 
emission properties. All have cool far-infrared colors. The optical classifications as LINERs 
are based on Veilleux et al. (1995) and Veilleux, Kim, & Sanders (1999). Table 2 summarizes 
the observing log. 

The 3.1-4.0 fim slit spectrum of UGC 5101 has already been presented in Imanishi et 
al. (2001). Descriptions of the observations and data analysis were offered in that Letter. To 
summarize, the observations were made at the IRTF on Mauna Kea, Hawaii, using the grism 
mode of NSFCAM (Shure et al. 1994). Although the wavelength coverage of the NSFCAM 
grism mode is 2.8-4.1 /im, results at A < 3.1 were not presented in Imanishi et al. (2001) 
because our main aim in that study was to investigate the 3.3 /im PAH emission and 3.4 /im 
carbonaceous dust absorption features, and because the Earth's atmospheric transmission 
at 2.8-3.1 /im is poorer than at A > 3.3 fim. However, after the publication of Imanishi 
et al. (2001), Spoon et al. (2002) reported the detection of 6.0 /im H 2 ice absorption in 
UGC 5101, and suggested that the 3.1 /im H2O ice absorption feature might be present in 
our 3.1-4.0 /xm spectrum. Therefore we carefully analyzed the spectrum at 2.8-3.1 fim to 
investigate the 3.1 fim H 2 ice absorption feature, whose profile is very broad, extending 
over at least A rest = 2.8-3.3 fim in the rest-frame (Smith, Sellgren, & Tokunaga 1989). 

The 2.8-4.1 fim spectra of the other ULIRGs were obtained with the Infrared Spectro- 
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graph and Camera (IRCS; Kobayashi et al. 2000) at the Subaru telescope on Mauna Kea, 
Hawaii. The sky was clear during the observations of these objects and the seeing at K was 
measured to be / /5-0 / /8 in full-width at half maximum. A O'.'9-wide slit and the L-grism were 
used with a 58-mas pixel scale. The achievable spectral resolution is ~140 at 3.5 /im. A 
standard telescope nodding technique, with a throw of 7 arcsec along the slit, was employed 
to subtract background emission. ULIRGs and corresponding standard stars were observed 
with an airmass difference of <0.1 to correct for the transmission of the Earth's atmosphere, 
and to provide flux calibration. 

Standard data analysis procedures were employed, using IRAF 4 . Initially, bad pixels 
were replaced with interpolated signals from the surrounding pixels. Bias was subsequently 
subtracted from the obtained frames and the frames were divided by a spectroscopic flat 
frame. Finally the spectra of ULIRGs and standard stars were extracted. Wavelength 
calibration was performed taking into account the wavelength-dependent transmission of 
the Earth's atmosphere. The spectra of ULIRGs were divided by the observed spectra of 
standard stars, multiplied by the spectra of blackbodies with temperatures appropriate to 
individual standard stars (Table 2), and then flux-calibrated. 

Appropriate binning of spectral elements was performed, particularly at Xobserved < 3.3 
jum in the observed frame, so as to give an adequate signal-to-noise ratio in each element. 
The resulting spectral resolution R at Xobserved < 3.3 fxm is ^100. Although the Earth's 
atmospheric transmission curve at 2.8-3.3 /im is highly wavelength-dependent if observed 
at high spectral resolution of R > 1000 (Figure la), it is fairly smooth at lower spectral 
resolutions, with R ^ 100 (Figure lb). Thus, even if the net positions of the target object 
and standard star on the slit differ slightly (on the sub-pixel scale) along the wavelength 
direction, the standard data analysis described above is expected to produce no significant 
spurious features in spectra with R ^ 100. 



3. Results 

Figure 2 shows flux-calibrated 2.8-4.1 /im spectra of the seven ULIRGs. Since our slit 
spectra cover physical scales larger than a few hundred pc in all the ULIRGs, emission from 
both AGNs and compact (less than a few hundred pc) starbursts should be fully detected. 



4 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc. (AURA), under cooperative agreement with the 
National Science Foundation. 
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All the ULIRGs clearly show the 3.3 /tm PAH emission. To estimate its flux, luminosity, 
and rest-frame equivalent width, we make the reasonable assumption that the profiles of the 
3.3 /tin PAH emission in these ULIRGs are similar to those of Galactic star-forming regions 
and nearby starburst galaxies; the main emission profile then extends between A resi = 3.24- 
3.35 /im (Tokunaga et al. 1991; Imanishi & Dudley 2000). The flux excess above the solid 
lines in Figure 2 should thus be ascribed to the 3.3 /itn PAH emission. Table 3 summarizes 
the results for the 3.3 /tin PAH emission in these seven ULIRGs. 

The spectrum of UGC 5101 shows clear 3.4 /im carbonaceous dust absorption with an 
observed optical depth of T3.4 = 0.65 (Imanishi et al. 2001). In other ULIRGs, there is no 
clear indication of this feature, partly because this absorption feature is intrinsically not so 
strong (r 3 . 4 /Av ~ 0.004-0.007; Pendleton et al. 1994). 

In addition to the 3.3 /tm PAH emission feature, there is another important common 
feature of the spectra of five ULIRGs (UGC 5101, IRAS 00188-0856, IRAS 03250+1606, 
IRAS 16487+5447, and IRAS 21329-2346): the continuum emission is concave in these 
spectra. At the shorter wavelength side of the 3.3 /im PAH emission, the continuum flux 
level initially decreases with decreasing wavelength, but then begins to increase again at 
A rest ~ 3.05 /im. In IRAS 00188-0856 and IRAS 03250+1606, the spectra become flat at 
the shortest wavelength parts (\ re .st < 2.7 /im). This behavior is quite similar to that of 
the infrared-luminous galaxy NGC 4945, in which 3.1 /tm H 2 ice absorption was clearly 
detected in an ISO spectrum (Spoon et al. 2000). 

Earlier ground-based L-band spectroscopy has revealed that many bright Galactic sources 
behind substantial columns of dust and molecular gas show broad 3.1 /tm H2O ice absorption, 
with a main profile peaking at A res t ~ 3.05 /tm and extending at least to \ re .st = 2.8-3.3 /tm 
(Smith et al. 1989; Smith, Sellgren, & Tokunaga 1993; Smith 1993). A recent IS O spectrum 
of NGC 4945 suggests that the absorption may extend to wavelengths as short as X re st ~ 2.5 
/tm (Spoon et al. 2000). In addition to this main profile, a weak longer- wavelength wing at 
\ rest = 3.3-3.7 /tm can also be distinguished (Smith et al. 1989, 1993; Smith 1993). The 
concave-shaped continuum emission found in the spectra of these five ULIRGs is consistent 
with the presence of this broad 3.1 /tm H 2 ice absorption. 

Some Seyfert 2 galaxies with lower infrared luminosities were observed on the same 
nights as these ULIRGs, but no signatures of strong 3.1 /tm H 2 ice absorption were found. 
The spectra of standard stars observed before these ULIRGs on the same nights were divided 
by other standard stars observed after, but the resulting spectra also never showed such 
signatures. The strong 3.1 /tm H 2 ice absorption signatures found only in the spectra of 
the ULIRGs are thus believed to be genuine and intrinsic to the ULIRGs, and not an artifact 
caused by insufficient cancellation of the Earth's atmospheric transmission. In fact, for UGC 
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5101 and IRAS 00188—0856, the clear detection of H 2 ice absorption at 6.0 /iin was also 
reported by Spoon et al. (2002). 

To estimate conservative lower limits for the optical depths of the broad 3.1 fim H 2 ice 
absorption feature, the concave quadratic dashed lines in Figure 2 are adopted as the lowest 
plausible continuum levels. In IRAS 23234+0946 and IRAS 23327+2913, the presence of 
strong 3.1 fim H 2 ice absorption is not clear in the case of the adopted concave continuum. 
For the remaining five ULIRGs, the observed optical depths of the 3.1 /im H 2 absorption 
(T3.1) as well as the 3.4 fim carbonaceous dust absorption (V3.4) are summarized in Table 4. 
These T3.1 values should be regarded as strict lower limits because: (1) linear continuum 
levels that connect data points at the shortest and longest wavelength parts of the spectra 
would give rise to larger values of r 31 , and (2) the data points at the shortest wavelength 
parts, used to determine continuum levels, may still be affected by the 3.1 fim H 2 ice 
absorption feature in some cases. Compared to UGC 5101, IRAS 00188—0856 and IRAS 
21329—2346, the presence of the broad H 2 ice absorption feature is apparently less clear 
in IRAS 03250+1606 and IRAS 16487+5447. However, as long as the spectral shapes of the 
3.3 fim PAH emission in these ULIRGs do not differ significantly from those observed so far 
in starburst galaxies, the continuum flux levels on the shorter wavelength side of the 3.3 fim 
PAH emission are clearly lower than on the longer wavelength side, which is consistent with 
the presence of the broad 3.1 fim H 2 ice absorption. The very broad emission-like features 
that extend between \ bserved = 3.6-3.9 fim and 3.4-3.8 fim in IRAS 03250+1606 and IRAS 
16487+5447, respectively, cannot be ascribed solely to the 3.3 fim PAH emission; instead, 
they are caused by a combination of 3.1 fim H 2 ice absorption, 3.3 fim PAH emission, and 
(possibly) 3.4 fim carbonaceous dust absorption. 



4. Discussion 

4.1. Energetic importance of weakly obscured starbursts 

The observed 3.3 fim PAH to infrared luminosity ratios for these seven ULIRGs are 
L3.3PAH/L1R ~ 1-3 xlO" 4 (Table 3), which are 3 to 10 times smaller than the ratios for 
starburst galaxies (L 3 . 3 pah/Lir ~ 10~ 3 ; Mouri et al. 1900; Imanishi 2002). It can therefore 
be concluded that weakly obscured (Ay < 20 mag) nuclear starbursts are not contributing 
significantly to the infrared luminosities of these ULIRGs. 

Since the L 3 3 pah values are based on slit spectroscopy measurements, the luminosities 
of extended (kpc scale) weakly obscured starbursts are not taken into account. However, as 
mentioned in §1, such extended starbursts have been shown to be energetically insignificant 
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in ULIRGs (Surace & Sanders 1999; Surace, Sanders, & Evans 2000; Fischer 2000). For UGC 
5101, Spoon et al. (2002) estimated, based on the ISO's large aperture measurement, a 6.2 
fxm PAH emission flux of 2.0 x 10~ 19 W cm -2 , so that its luminosity is ^6.2PAh ~ 6.5 x 10 42 
ergs s" 1 in our adopted cosmology. The resulting L6.2PAh/-^ir ~ 1-6 x 10~ 3 , which is only 
~ 1/4 of the ratio typically seen in starburst galaxies (6 x 10~ 3 : Fischer 2000). Thus, weakly 
obscured extended starbursts are unlikely to be the dominant energy source in UGC 5101, 
either. The dominant energy sources should be compact and heavily obscured; for the five 
ULIRGs, the presence of such buried energy sources is supported by the detection of strong 
absorption features in our 2.8-4.1 /im spectra. 



4.2. Geometrical requirements on the buried energy sources 

The detection of strong (Y3.1 > 0.6) 3.1 /im H 2 ice absorption can, in principle, pro- 
vide stringent constraints on the nature of the buried energy sources at the nuclei of the 
five ULIRGs. The 3.4 /im carbonaceous and 9.7 /im silicate dust absorption features are 
detectable if an adequate amount of dust in the diffuse inter-stellar medium is present in 
front of the continuum-emitting energy sources. However, in order for H 2 ice absorption 
to be detected, dust grains covered with an ice mantle must be present in the foreground. 
Such a situation can arise only when the dust grains are located in molecular clouds, and are 
sufficiently shielded from the UV radiation from any energy sources either inside or outside 
the molecular clouds (Whittet et al. 1988). In other words, the detection of strong H 2 ice 
absorption requires that a large amount of dust in molecular clouds be present in front of 
the continuum-emitting energy sources. Foreground molecular clouds in the host galaxies 
could be responsible for the strong 3.1 /im H 2 ice absorption detected toward the nuclei 
if the host galaxies were seen from nearly edge-on. However, there is no strong evidence 
that the host galaxies of these ULIRGs are edge-on (Murphy et al. 1996; Scoville et al. 2000; 
Kim et al. 2002). Spoon et al. (2002) reported that 6.0 /im H 2 ice absorption is systemat- 
ically stronger in ULIRGs than in less infrared-luminous Seyfert 2 galaxies. These Seyfert 2 
galaxies also generally show no detectable or weak (r 3 1 < 0.2) 3.1 /im H 2 ice absorption 
(Imanishi et al. 2002, in preparation). The systematically larger optical depths of H 2 ice 
absorption toward the nuclei of ULIRGs can reasonably be explained by a larger amount 
of highly concentrated molecular gas (and dust) in their nuclei (Sanders & Mirabel 1996). 
We thus regard it as likely that the strong (r 31 > 0.6) 3.1 /im H 2 ice absorption detected 
toward the nuclei of these five ULIRGs are caused by nuclear concentrations of molecular 
gas and dust. Note that, although the host galaxy of NGC 4945 is viewed edge-on (Abies et 
al. 1987), a significant fraction of the strong 3.1 /im H 2 ice absorption (T3.1 ~ 0.9; Spoon 
et al. 2000) is likely to originate in a highly concentrated nuclear component of ice-covered 
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dust grains, on a scale of less than a few hundred pc (Brock et al. 1988). 

In the case of normal starbursts, it seems reasonable to assume that the continuum- 
emitting sources and molecular gas and dust are spatially well mixed (Soifer et al. 2000). 
A detailed infrared study of the nearby starburst galaxy M82 suggests that the mixed-dust 
geometry is a better fit to the observed data than a foreground screen dust distribution 
(McLeod et al. 1993). Since the oscillator strengths of species are fixed, there are upper 
limits on the absorption optical depths that can be produced by the mixed-dust geometry 
(unless unusual abundances are considered). This is because the foreground, less-obscured, 
and therefore less-attenuated, emission (which shows only weak absorption features) will 
dominate the observed fluxes. In practice, nearby starburst galaxies such as M82 and NGC 
253 show weak (r ~ 0.2) or undetectable absorption by H 2 ice (Sturm et al. 2000) or dust 
grains (Roche et al. 1991; Dudley & Wynn- Williams 1999; Imanishi & Dudley 2000). This 
issue will be discussed quantitatively below. 

In Galactic molecular clouds, the optical depth of the H 2 ice absorption and dust 
extinction (Ay) are correlated, with a relation given by r 31 = a x (Ay — Ay ). The value 
of a has been found to be almost constant, ~0.06, in Galactic molecular clouds, which 
exhibit various different levels of star-forming activity (Smith et al. 1993; Tanaka et al. 1990; 
Murakawa, Tamura, & Nagata 2000). This is because it reflects the abundance of H2O 
ice relative to dust grains in the regions where an ice mantle can survive; it is therefore 
unlikely to differ substantially from the Galactic value, even in starburst galaxies. The 
Ayo value in the above relation is the threshold required for dust grains to be sufficiently 
shielded from the ambient UV radiation. This value can vary dramatically, from 2-6 mag in 
quiescent molecular clouds (Taurus and Serpens; Whittet et al. 1988; Eiroa & Hodapp 1989; 
Murakawa et al. 2000) to 15 mag in active star-forming molecular clouds (p-Oph; Tanaka 
et al. 1990). In starburst galaxies with a more intense UV radiation field than the p-Oph 
molecular clouds, Ay is almost certainly larger than 15 mag, but its exact value is unknown. 
Thus the relation 



where / is the fraction of dust that is covered with an ice mantle. 

Figures 3a and 3b show schematic diagrams of a mixed-dust geometry. In this geometry, 
the observed flux I(t„) is given by 



T3.1 = 0.06 x (A v - A yo ), 



(1) 



is adopted in starburst galaxies, where Ayo > 15 mag, or 



t 31 = 0.06 x / x A 



(2) 




= h x 



1 -e- 



(3) 



T, 
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where Jo is an unattenuated intrinsic flux and r v is the optical depth at each wavelength, 
which takes different values inside and outside the absorption features. For 3-4 /xm con- 
tinuum emission outside the 3.1 /im H 2 ice and 3.4 /im carbonaceous dust absorption 
features, 

T cont = 0.06 xA v (4) 

(Rieke & Lebofsky 1985; Lutz et al. 1996), where, by definition, Ay = 1.08 x ry. At the 
peak wavelength of the 3.1 /jm H 2 ice absorption feature, 

T ice = 0.06 X Ay + 0.06 X f X Ay, (5) 

where the first and second terms give the flux attenuation for 3-4 /iin continuum and the 
H 2 ice absorption feature, respectively. Thus, the optical depth of the H 2 ice absorption 
feature in an observed 3-4 /im flux is 

-3.i = ln[— x— T ^-\ (6) 

1 _ p-0.06x.4v 

= M(l + /) 1 _ e - .o,xAvx(l + /) ] ( ? ) 

This last formula implies that, for a fixed /-value, r 3 i increases with increasing Ay, but 
saturates at a certain value. The /-value can be constrained based on observations of the 
nearby starburst galaxy M82. Based on a mixed-dust model, the dust extinction toward the 
farthest parts of the starburst region is estimated to be Ay ~ 55 mag (McLeod et al. 1993). 
The observed optical depth of the 3.1 /im H 2 ice absorption is t 31 ~ 0.2 (Sturm et al. 
2000). Based on these values, / ~ 0.3 is obtained. In other words, in the starburst galaxy 
M82, ~ 30% of dust grains are covered with an ice mantle. 

If the cores of ULIRGs are simply a scaled-up version of normal starbursts with a 
mixed-dust geometry and the /-value is ~0.3, then r 31 should be ~0.2. However, the /- 
value is strongly dependent on the fraction of active star-forming regions relative to quiescent 
molecular gas in molecular clouds. If the cores of ULIRGs are powered by starbursts with 
a mixed-dust geometry, the /-value in ULIRGs should be even smaller than in less infrared 
luminous starburst galaxies (Soifer et al. 2000). Here, however, / ~ 0.3 is adopted for 
ULIRGs in order to determine a conservative upper limit for r 3 .i within the framework of a 
mixed-dust geometry. Even if this conservative value is adopted, Equation (7) shows that 
T3.1 cannot be larger than 0.3. Therefore, it can be concluded that the observed r 31 values of 
>0.6 in these five ULIRGs are clearly incompatible with mixed-dust geometry and f ^ 0.3. 
When the contribution from any weakly obscured starburst emission to the observed 3-4 
/im fluxes is subtracted, the r 3 .i values toward the buried energy sources will be even larger, 
making this conclusion more robust. 
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The same argument is applicable to the 3.4 /j,m carbonaceous dust absorption in UGC 
5101 (Imanishi et al. 2001). For the 3.4 /xm carbonaceous dust absorption feature, 

Tcarbo = 0.06 X Ay + (0.004 ~ 0.007) X (1 - /) X Ay (8) 

is obtained (Pendleton et al. 1994), where the term (1 — /), the fraction of dust grains 
not covered by an ice mantle, is included because the 3.4 /zm carbonaceous dust absorption 
feature may be suppressed if the dust grains are covered with an ice mantle (Mennella et al. 
2001). Based on the same arguments as above, it can be quantitatively demonstrated that 
r 3 4 cannot be larger than 0.2 in the mixed-dust geometry. The observed large optical depth 
of T3.4 > 0.6 in UGC 5101 (Imanishi et al. 2001) is again incompatible with a mixed-dust 
source geometry. 

In order to explain the strong absorption features in these five ULIRGs, dust in a 
foreground screen must be obscuring the 3-4 /zm continuum-emitting energy sources. This 
configuration requires that the energy sources be centrally concentrated and deeply buried 
in nuclear dust and molecular gas, as shown in Figure 3c. 

4.3. What are the centrally-concentrated energy sources? 

4.3.1. The case for UGC 5101 

The centrally-concentrated nature of the buried energy source in UGC 5101 is naturally 
explained by a buried AGN (Soifer et al. 2000). Furthermore, as mentioned in §1, the 
observed rest-frame equivalent widths of the 3.3 /im PAH emission in starburst galaxies 
should be little changed by dust extinction. The significantly lower equivalent width in 
UGC 5101 (Table 3) than in starbursts (~120 nm; Moorwood 1986; Imanishi & Dudley 
2000) also supports the presence of AGNs that produce strong 3-4 /iin continuum emission, 
but essentially no PAH emission. 

When our data are combined with ground-based high spatial resolution mid-infrared 
(8-25 /im) imaging data (Soifer et al. 2000), we can give more insight into the nature of the 
buried energy source. Soifer et al. (2000) estimated the size of the spatially-unresolved core 
of UGC 5101 at 12.5 /im to be <200 pc. If the attenuation of the mid-infrared flux from 
the spatially-unresolved compact core in UGC 5101 is negligible, the surface brightness of 
the core emission becomes >1.0 x 10 12 L Q kpc~ 2 (Soifer et al. 2000). However, the detection 
of strong 3.4 /jm carbonaceous dust and 3.1 fim H 2 ice absorption toward the nucleus of 
UGC 5101 requires that 

1. the actual size of the buried centrally-concentrated energy source is substantially 
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smaller than that of the dust and molecular gas, and 

2. the dust extinction toward the buried energy source is Ay > 115 mag (Imanishi et al. 
2001). 

On point (1), since dust around a centrally-concentrated energy source shows a temper- 
ature gradient (see §1), the size of the cooler dust region, which is also responsible for the 
foreground dust extinction, is larger than that of the ~250 K dust traced by the 12.5 fim 
emission (<200 pc). However, the size of the cool dust emission (traced by the 1.3 mm emis- 
sion) at the core of the ULIRG Arp 220 has been shown to be compact (<200 pc; Sakamoto 
et al. 1999). Furthermore, the sizes of the high density molecular gas regions at the cores of 
ULIRGs have been estimated to be no more than ~500 pc in diameter (Solomon et al. 1992). 
Thus, it seems reasonable to assume that the actual size of the buried centrally-concentrated 
energy source in UGC 5101 is only a fraction of this, ~100 pc or less. On point (2), the 
dust extinction of Ay > 115 mag means that the flux attenuation at 12.5 /xm is more than 
factor of 25 (Rieke & Lebofsky 1985). When these corrections are taken into account, the 
actual surface brightness of the centrally-concentrated buried energy source in UGC 5101 
is > 1O 14 L kpc~ 2 . This exceeds the peak surface brightness of star clusters in starburst 
galaxies (~5 x 1O 13 L kpc -2 ; Meurer et al. 1997; Soifer et al. 2000) by a factor of more 
than two. The high surface brightness is best explained by the presence of an energetically 
important buried AGN in UGC 5101. 

If a buried AGN is present at the nucleus of UGC 5101, the dust is expected to show 
a strong temperature gradient (see §1). Are the available data consistent with its presence? 
Soifer et al. (2000) could not explicitly find the signature of a temperature gradient based 
on their imaging data of UGC 5101. Here, it is searched for using spectroscopic data. As 
mentioned in §1 and Imanishi (2000), if a strong temperature gradient is present in the dust, 
the optical depth of the 3.4 /iin carbonaceous dust absorption feature probes deeper into the 
dust distribution than that of the 9.7 silicate dust absorption. Thus, the T3.4/V9.7 ratio 
toward the nucleus of UGC 5101 is expected to be significantly larger than the value seen 
when no significant temperature gradient is present (r 3A /rg, 7 ~ 0.06-0.07; Roche & Aitken 
1984, 1985; Pendleton et al. 1994). 

After subtracting the contribution from the weakly obscured starburst emission to the 
observed 3-4 /iin flux, the T3.4 value toward the buried AGN at the core of UGC 5101 was 
estimated to be ~0.8 (Imanishi et al. 2001), while only a lower limit for the observed r 9 7 
value (>1.5) was given by Spoon et al. (2002) based on a measurement with a large aperture 
using ISO. Although the r 9 . 7 value is a lower limit, its actual value is unlikely to exceed 2.0, 
judging from the actual spectrum of UGC 5101 (Spoon et al. 2002; their Figure 8). Thus 
the T3. 4 /t 9 . 7 ratio is > 0.4. The actual r 9 . 7 value toward the buried AGN in UGC 5101 could 
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be larger than the observed value, because weakly obscured starburst emission is likely to 
contribute significantly to the observed ~10 /zm ISO flux. However, since the equivalent 
width of the 6.2 /zm PAH emission in the ISO spectrum of UGC 5101 is less than half that of 
starburst galaxies (Spoon et al. 2002), at least half of the observed continuum flux at >5 /xm 
should originate in the buried AGN, implying that the T3.4/V9.7 ratio toward the buried AGN 
is > 0.2. Thus, these spectroscopic data are consistent with the presence of a temperature 
gradient of dust at the core of UGC 5101. 

If the dust around an AGN has a torus-like geometry, low-density clumpy clouds which 
are directly illuminated by the central AGN can exist at relatively large distances from the 
AGN along the axis of the torus-geometry dust (the so-called the narrow line regions or 
NLRs). NLRs produce the strong high-excitation forbidden emission lines seen in AGN. 
However, in a buried AGN, X-ray dissociation regions (Maloney et al. 1996), rather than 
NLRs, are expected to be produced around the AGN, so that high-excitation forbidden 
emission lines will be weak, consistent with the observed data of UGC 5101 (Genzel et al. 
1998). 

Solomon et al. (1992) and Solomon (2001) found that the infrared to HCN (J=l-0) 
luminosity ratios in ULIRGs, including UGC 5101 (Gao & Solomon 2002), roughly follow 
the relation established in less infrared-luminous starburst galaxies. By regarding the HCN 
( J=l-0) luminosities as a tracer of the mass of high density (>10 4 cm -3 ) molecular gas where 
star- formation actually occurs, Solomon et al. (1992) argued that the infrared luminosities 
of ULIRGs are powered by starbursts. However, in the luminous inner regions (central few 
hundred pc) of ULIRGs, the gas pressures - and therefore molecular gas densities - are 
inevitably extremely high, whether the energy sources are starbursts or AGNs (Barvainis et 
al. 1997; Maloney 1999). The relative constancy of the infrared to HCN luminosity ratio 
from less infrared luminous starburst galaxies to ULIRGs simply indicates that much of the 
infrared emission arises from dust in dense molecular clouds, which is inevitable whatever 
the sources of the luminosities of ULIRGs. 



4.3.2. The other ULIRGs 

For the two ULIRGs which do not show strong 3.1 /im H 2 ice absorption (IRAS 
23234+0946 and IRAS 23327+2913), there is no strong observational evidence for the pres- 
ence of buried AGNs. However, for the other four ULIRGs which show strong 3.1 /im H2O ab- 
sorption (IRAS 00188-0856, IRAS 03250+1606, IRAS 16487+5447, and IRAS 21329-2346), 
the centrally-concentrated nature of their buried energy sources and systematically (roughly 
a factor of two) lower-equivalent-width PAH emission than starbursts (Table 3) are also nat- 
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urally explained by the presence of energetically important, buried AGNs. Assuming that 
approximately half of the observed 3-4 /im fluxes originate in buried AGNs, the r 3 i and r 3 4 
toward the buried AGNs are twice as large as the observed values. These r 3 i and r 3 . 4 values 
reflect the column density of foreground dust, toward the buried AGNs, covered with and 
without an ice mantle, respectively (§4.2). Based on the relation of 73.1 /Ay ~ 0.06, r 3 .4/Ay 
~ 0.004-0.007 (§ 4.2), and A 3 _^ m / Ay ~ 0.06 (§ 1), and on a foreground screen dust model 
(§ 1), it is found that the dereddened 3-4 //m dust emission luminosities powered by the 
AGNs could be comparable to the infrared luminosities of these four ULIRGs. 

It is sometimes argued that ULIRGs with cool far-infrared colors cannot be powered by 
AGNs because AGNs show warm far- infrared colors (e.g., Downes & Solomon 1998). Indeed, 
many known Seyfert 2 galaxies show warm far-infrared colors (de Grijp, Miley, & Lub 1985). 
However, as mentioned in §4.2, the systematically higher detection rate of H 2 ice absorption 
in ULIRGs, as compared to less infrared luminous Seyfert 2 galaxies, is naturally explained 
by a higher column density of nuclear dust in ULIRGs. This explanation is supported by 
imaging data at ~10 /xm: while emission from the nuclei of less infrared luminous Seyfert 2 
galaxies is argued to be optically thin at ~10 jum (Alonso-Herrero et al. 2001), that from the 
nuclei of ULIRGs is estimated to be optically thick (Soifer et al. 2002). The higher column 
densities of dust around buried AGNs in ULIRGs can naturally produce infrared spectral 
energy distribution with cooler color than in less infrared luminous Seyfert 2 galaxies. It is 
quite possible that many other ULIRGs with cool far-infrared colors contain energetically 
important AGNs deeply buried in dust. 

Takeuchi, Yoshikawa, & Ishii (2002) have found that a 60 /zm luminosity function in 
the local universe, derived based on ~15000 galaxies, shows a bump at L 60 ^ m ~ 10 12 L o , and 
suggested that AGN activity begins to be energetically important at L 6OMm ^lO 12 L . This 
result also favors the buried AGN scenario in the majority of nearby ULIRGs. 

5. Summary 

Strong (t 3 .i > 0.6) 3.1 /im H 2 ice absorption has been detected toward the nuclei of 
five ULIRGs: UGC 5101, IRAS 00188-0856, IRAS 03250+1606, IRAS 16487+5447, and 
IRAS 21329—2346. All these ULIRGs have cool far-infrared colors and are classified opti- 
cally as LINERs. This detection requires a foreground screen dust geometry, which is best 
explained if centrally-concentrated, compact energy sources are buried in nuclear dust and 
molecular gas. For all the ULIRGs, the observed equivalent widths of the 3.3 /iin PAH emis- 
sion are significantly lower than those of starburst galaxies. These observational data are 
naturally explained by the presence of AGNs deeply buried in dust, and support our previ- 



-16- 



ous arguments that LINER ULIRGs with cool far-infrared colors may possess energetically 
important, buried AGNs. 
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Table 1. Observed LINER ULIRGs and their far-infrared emission properties. 



Object 


Redshift 


fi2 a 


f 25 a 


r a 
160 


boo 


log Lfir b 


log L IR c 


f25/f60 


d 






(Jy) 


(Jy) 


(Jy) 


(Jy) 


(ergs s _1 ) 


(ergs s _1 ) 






UGC 5101 


0.040 


0.25 


1.03 


11.54 


20.23 


45.5 


45.5 


0.09 (cool) 


IRAS 00188-0856 


0.128 


<0.12 


0.37 


2.59 


3.40 


45.8 


45.9-46.0 


0.14 (cool) 


IRAS 03250+1606 


0.129 


<0.10 


<0.16 


1.38 


1.77 


45.6 


45.6-45.7 


<0.12 (cool) 


IRAS 16487+5447 


0.104 


<0.08 


0.20 


2.88 


3.07 


45.6 


45.7 


0.07 (cool) 


IRAS 21329-2346 


0.125 


<0.08 


<0.16 


1.65 


2.22 


45.6 


45.6-45.7 


<0.10 (cool) 


IRAS 23234+0946 


0.128 


<0.06 


<0.20 


1.56 


2.11 


45.6 


45.6-45.7 


<0.13 (cool) 


IRAS 23327+2913 


0.107 


<0.06 


0.22 


2.10 


2.81 


45.6 


45.6-45.7 


0.10 (cool) 



a fi2) f25; f60) an d boo are IRAS FSC fluxes at 12 /an, 25 ^m, 60 ^m, and 100 /um, respectively. 

b Logarithm of far-infrared (40-500 fim) luminosity in ergs s _1 calculated with Lfir = 1.4 x 1.5 x 10 39 x 
D(Mpc) 2 x(2.58 x / 60 + /ioo) ergs s" 1 (Sanders & Mirabel 1996). 

c Logarithm of infrared (8—1000 fim) luminosity in ergs s" - 1 calculated with L m = 2.1 x 10 39 x D(Mpc) 2 
x (13.48 x /12 + 5.16 x / 25 + 2.58 x / 60 + / 100 ) ergs s" 1 (Sanders & Mirabel 1996). 

d f 2 5/f60 < (>) 0.2 are called cool (warm) (Sanders et al. 1988). 
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Table 2. Observing log 



Object 


Date 


Integration 


P.A. a 




Standard Stars 






(UT) 


(Min) 


(°) 


Name 


L-mag 


Type 


T cff (K) 


UGC 5101 


2001 April 9 


18 


90 


HR 4112 


3.4 


F8V 


6200 


IRAS 00188-0856 


2002 August 19 


25 


90 b 


HR 72 


5.0 


GOV 


5930 


IRAS 03250+1606 


2002 August 19 


23.3 





HR 763 


4.3 


F7V 


6240 


IRAS 16487+5447 


2002 March 28 


36 


70 c 


HR 5949 


6.3 


A0V 


9480 


IRAS 21329-2346 


2002 October 24 


32 


90 


HR 7898 


4.5 


G8V 


5400 


IRAS 23234+0946 


2002 October 24 


32 


120 d 


HR 8653 


4.6 


G8IV 


5400 


IRAS 23327+2913 c 


2002 October 24 


24 


90 


HR 8955 


5.1 


F6V 


6400 



a Position angle of the slit. 0° corresponds to the north-south direction. Position angle increases with 
counter-clockwise on the sky plane. 

b To avoid the contamination from the emission of the southern foreground star (Kim, Veilleux, & 
Sanders 2002). 

c To attempt to detect signals from both components of the double nucleus of this object, which are 
separated by ~3 arcsec in the X-band (2.2 ^m) image (Murphy et al. 1996). Only the western primary 
nucleus (Murphy et al. 2001) was clearly detected in our 2.8-4.1 /j,m spectrum. 

d To attempt to detect signals from both components of the double nucleus of this object, which are 
separated by ~5 arcsec in the K-band image (Kim et al. 2002). Only the western nucleus was clearly 
detected. 

c The southern primary nucleus (Trung et al. 2001) was observed. 
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Table 3. The properties of the 3.3 fim PAH emission. 



(1) 


(2) 


(3) 


(4) 


(5) 


Object 


f3.3PAH 


L3.3PAH 


L3.3PAH/L1R 


rest EW3.3PAH 




(xl0~ 14 ergs s _1 cm~ 2 ) 


(xl0 41 ergs s- 1 ) 


(xl0~ 4 ) 


(nm) 


UGC 5101 


16+3 


5.2+1.0 


~1 


~35 


IRAS 00188-0856 


3.6+0.4 


13+1 


~1 


~50 


IRAS 03250+1606 


3.2+0.4 


12+1 


~3 


~80 


IRAS 16487+5447 


2.6+0.4 


6.2+1.0 


~1 


~70 


IRAS 21329-2346 


1.2+0.2 


4.2+0.8 


~1 


~50 


IRAS 23234+0946 


2.6+0.3 


9.9+1.0 


~2 


~70 


IRAS 23327+2913 


1.9+0.3 


4.8+0.8 


~1 


~50 


Note. — Col. (1): 


Object name. Col. (2): Observed flux of the 3.3 (im PAH emission. Col. (3): 


Observed luminosity of the 3.3 fim. PAH emission. Col. (4): Observed 3.3 /im PAH-to-infrared 


luminosity ratio. Col. (5): Rest-frame equivalent width of the 3.3 /j,m PAH emission. Those for 




starbursts are 


~120 nm (see text). 




Table 4. The observed optical depths of the 3.1 H 2 


ice and 3.4 jum 


carbonaceous 


dust absorption 


features for the five ULIRGs with strong 


3.1 jum H 2 ice 


absorption. 




Object ol 


Dserved T3.1 observed r 3 4 






UGC 5101 


>0.8 


0.65 






IRAS 00188-0856 


>1.6 


<0.2 






IRAS 03250+1606 


>0.6 


<0.1 






IRAS 16487+5447 


>0.8 


<0.35 






IRAS 21329-2346 


>0.9 


<0.2 





AA 



0.002/xm 



AA 



0.064/xm 



o 



o 



3 



3.5 

wavelength (jim) 



4 



3 



3.5 



4 



wavelength (/Aim) 



Fig. 1. — The Earth's atmospheric transmission at 2.8-4.1 /im at the summit of 
Mauna Kea under good weather conditions. Adopted from the UKIRT webpage 
at http:/ /www .jach.hawaii.edu/JACpublic/UKIRT/astronomy/calib/indexl-tx.html. (a): 
Measured with a spectral resolution of 1500-2000 at 3-4 /im (AA = 0.002 /im). (b): Mea- 
sured with a spectral resolution of 50-60 at 3-4 /im (AA = 0.064 /im), roughly the same as 
the spectral resolution of the seven ULIRGs at Xobserved < 3.3 /im. Note that the transmission 
becomes worse with increasing precipitable water, particularly at 2.8-3.3 /im. 
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Fig. 2. — 2.8-4.1 /im spectra of the seven ULIRGs. The abscissa and ordinate are the 
observed wavelength and the flux in F\, respectively. The dashed concave lines are the 
continuum levels used to estimate conservative lower limits for the optical depths of the 3.1 
/im H 2 ice absorption. The solid lines are the continuum levels used to estimate the values 
of the 3.3 /im PAH emission and the 3.4 /im carbonaceous dust absorption. 
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Fig. 3. — Possible geometries of the energy sources, molecular gas, and dust. In all cases, the 
dotted areas indicate the regions where an ice mantle can survive (that is, dust is covered with 
an ice mantle). The star symbols represent starbursts. The stars outside the outermost circles 
represent weakly obscured starbursts; only the side near the observer is shown, (a): Starburst 
activity is homogeneously mixed with dust and molecular gas. The fraction of active star- 
forming regions relative to quiescent molecular gas is modest, (b): The same geometry as 
(a), but active star-forming regions occupy the majority of the volume of molecular gas. 
(c): The buried energy sources are more centrally concentrated than the dust and molecular 
gas. The black dot represents an AGN. This is the most plausible geometry to account for 
the observed 2.8-4.1 /im spectral shapes of the five ULIRGs with strong 3.1 /jm H 2 ice 
absorption (see the text). 



